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Abstract 

We report first principles linear response calculations on NaCoC>2- Phonon frequencies and 
eigenvectors are obtained throughout the Brillouin zone for two geometries with different Na site 
occupancies. While most of the phonon modes are found to be unsensitive to the Na site occupancy, 
there are two modes dominated by out-of-plane vibrations of Na giving very different frequencies for 
different geometries. One of these two modes, the mode, is infrared-active, and can be used as a 
suitable sensor of Na distribution /ordering. The longitudinal-transverse splitting of the zone-center 
optical-mode frequencies, Born effective charges and the dielectric constants are also reported, 
showing considerable anisotropy. The calculated frequencies of Raman-active modes generally agree 
with the experimental values of corresponding Na de-intercalated and/or hydrated compounds, 
while it requires better experimental data to clarify the infrared-active mode frequencies. 

PACS numbers: 63.20.Dj, 71.20.Ps, 74.25.Kc 
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I. INTRODUCTION 



Since the discovery of the superconductivity in Na^CoCVyH^O,— significant interest 
has been focused on this novel material, because there are some intriguing similarities 
between this cobalt oxide superconductor and copper oxide superconductors.— A num- 
ber of evidences in experiments suggest that the superconductivity in this material is 
unconventional.—^ Theoretical models such as resonating valence bond (RVB ) 6 ' 7 ' 8 ' 9 and 
spin-triplet superconductivit y 10 ' 11 ' 12 ' 13 with strong magnetic quantum fluctuations are also 
proposed. Although now it is a consensus that the cobalt oxide superconductor is not simply 
a BCS superconductor, some experiments still suggest the possible strong lattice coupling 
to the electronic density of states in this compound . 14 ' 15 It is thus important to fully char- 
acterize the lattice dynamical properties of Na x Co02-yH20 and related materials. 

While there are some Raman , 16 ' 17 infrared , 14 ' 18 and neutron scattering^ experiments 
involving the lattice dynamics of the cobalt oxide compound, there is no first principles 
theoretical report found in the literature. In this article, we report a density functional 
perturbation theory (DFPT) 20 ' 21 calculation on NaCoC>2, the intrinsic (without hydration 
or de-intercalation of Na) insulating phase of the cobalt oxide superconductor. As shown 
in Figures la and lb, NaCoC>2 has a hexagonal structure (space group #194, P63/HIHIC) 
consisting of C0O2 layers of edge sharing CoOg octahedra and Na layers with two partly 
occupied Na Wickoff sites. To accommodate for the partial occupation of the two Na sites, 
we consider two geometries of NaCoC>2 with one site fully occupied and the other site empty, 
namely geometry A (Figure la) with the 2b site occupied and geometry B (Figure lb) with 
the 2d site occupied. 



II. COMPUTATIONAL METHOD 



The present results have been obtained through the use of the ABINIT code,— a common 
project of the Universit Catholique de Louvain, Corning Incorporated, and other contribu- 
tors (URL |http:/7 www.abinit.org), based on pseudopotentials and plane waves. It relies on 



an efficient fast Fourier transform algorithm for the conversion of wave functions between 
real and reciprocal spaces, on the adaptation to a fixed potential of the band-by-band con- 
jugate gradient method^ 3 -, and on a potential-based conjugate-gradient algorithm for the 



2 



determination of the self-consistent potential.— Technical details on the computation of re- 
sponses to atomic displacements and homogeneous electric fields can be found in Ref. 0, 
and the subsequent computation of dynamical matrices, Born effective ch arg es, dielectric 
permittivity tensors, and interatomic force constants was described in Ref. 125 . 

Troullier-Martins norm conserving pseudopotentials 2 ^ are used, with Teter 
parametrization 27 of the Ceperley-Alder exchange-correlation potential. The kinetic 
energy cutoff of the plane wave bases is 50 hartree. A uniform 6x6x2 /c-point mesh is 
used in Brillouin-zone integrations, which gives well converged total energies and phonon 
frequencies. 



III. RESULTS AND DISCUSSION 

A. Geometries and Electronic Structures 



As a first step, we optimize the geometry and calculate the corresponding electronic 
structure. During the geometry optimization, the cell parameters are fixed to the experi- 
mental values of Na .74CoO 2 (a = 2.84A, c = 10.811A).— We notice that neutron scattering 
experiment o 28 i 29 i 30 for samples with different Na concentrations give similar cell parameters. 
Considering the space group symmetry, the only freedom to be relaxed is the Co-0 layer 
space d, which turns out to be 1.897 and 1.899 A for geometries A and B respectively. The 
electronic band structure and density of states (DOS) based on the optimized structure of 
geometry A are shown in Figure lc, a gap of 0.94 eV (from 0.60 to 1.54 eV) between the 
valence band and the conduction band can be clearly identified. Little difference on the 
electronic structure of the two geometries is found. This similarity strongly indicates that 
Na contributes to the electronic structure only by doping electrons to the CoO layers. 



B. Phonon Dispersions 

To map the phonon dispersion curves throughout the Brillouin zone, the dynamical ma- 
trices are obtained on a uniform 6x6x2 grid of q points, and real-space force constants 
are then found by Fourier transform of the dynamical matrices. The dynamical matrix at 
an arbitrary wave vector q can then be computed by an inverse Fourier transform. The 
acoustic sum rule is applied to force the three acoustic phonon frequencies at T equal to zero 
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strictly as being implied by the translation symmetry. The calculated phonon dispersions 
and corresponding DOS for geometries A and B are shown in Figure 2. We see first that 
both structures are stable. Second, we notice that the phonon modes are separated to two 
groups in frequency, namely the soft group with frequencies below than about 400 cm -1 
and the hard group with frequencies between 450 to 650 cm -1 . As a whole, the phonon 
dispersions for the two geometries with different Na site occupancies are similar, but there 
are also some minor differences exist for the two geometries in the phonon band structure, 
unlike the nearly identical electronic band structures. First, as we will discussed in detail 
below, there are some soft modes giving significant different frequencies near zone center for 
the two geometries. Secondly, we notice that the dispersion along V — A for geometry B is 
larger than that for geometry A. 

C. zone-center phonon modes 

The zone-center phonon modes are of special importance, since they can be obtained by 
various of experimental methods. In Table I, we list the frequencies, symmetries and vi- 
bration modes of the zone-center phonons. The triply-degenerated acoustic mode with zero 
frequency is not listed. The frequency differences between the two geometries are generally 
small, especially for the Raman-active phonon modes. However there are an infrared- active 
A 2u mode and a silent B ls mode giving very different frequencies for geometry A and ge- 
ometry B, and both of these two modes are mainly contributed by the Na out-of-plane 
vibrations. The other A 2u and B l9 modes (modes 7 and 14) also involve the out-of-plane 
vibrations of Na, but we notice that the vibrations of Co and/or O in these two modes 
are much stronger than or at least comparable to the Na vibrations. Therefore they give 
relatively small frequency differences for different Na site occupancies. We find that the 
frequency difference of Raman-active A lg mode of O out-of-plane vibrations is rather small, 
in contrast to the conclusion of the shell model calculations by Lemmens et al..-^ They found 
a pronounced frequency dependence of the A ls mode, and argued that it would be used as 
a very susceptible sensor of Na distribution/ordering. However, in our first principles cal- 
culations, the infrared- active A 2u mode (node 6) in the soft group may be more suitable to 
act as the sensor mode. Another difference between T phonon modes of the two geometries 
is the frequency order of the first E 2s and B 2n modes (modes 1 and 10). 
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Since the systems we studied here are polar materials, we also get the longitudinal- 
transverse optical mode (LO-TO) splittings by considering response to the macroscopic 
electric field. As shown in Table I, the LO-TO splittings are found only in the infrared- 
active 2E lu + 2A 2u modes. The biggest LO-TO splitting is found in the hard E ltl mode, 
with Au = ujlo — ^to equal to 52.4 and 48.6 cm -1 for geometries A and B respectively. 
Contrastingly, the soft Ei u mode is split by only 5.1 and 5.9 cm _1 respectively. Both these 
two Ei u modes are "layer sliding" modes, and the hard one is mainly contributed by the 
Co layers sliding against the O layers. While for the soft E lu mode, the sliding of the Na 
layers is about an order stronger than the relative sliding of Co and O layers. Therefore our 
results follow the trend that displacements modulating the "covalent" bonding produce the 
largest LO-TO splitting, as being pointed out by Lee et al.— . 

We also get the Born (dynamical) effective charge tensor. In hexagonal symmetry, it is 
diagonal and reduces to two values Z xx = Z yy — Z// and Z zz = Z±. As listed in Table 
II, these charges show considerable anisotropy in NaCo02- We notice that Z± is much 
larger than Z// for Na, which is consistent with the fact that the LO-TO splitting of the 
soft in-plane E\ u mode is much smaller than that of the soft out-of-plane A 2u mode, with 
both modes being dominated by the vibrations of Na. Despite the anisotropy, we define the 
average effective charge Z as \TrZ, which gives about +1, +2 and -| for Na, Co, and O, 
respectively. The dynamical charges of Co and O are both smaller than their nominal ionic 
charges. In addition, the static electronic dielectric constants are calculated to be = 9.72, 

= 4.34 for geometry A and e{£ = 9.78, = 4.27 for geometry B respectively. 

There are some experimental data on the zone-center phonon frequencies of Na x Co02 
or Na a; Co02-yH20. In principle, these frequencies can not be directly compared with our 
results, since the systems we studied here is NaCo02- But as we have seen, the Na site 
occupancy may not affect the frequencies much, especially for the high-frequency in-plane 
vibration modes. So we also list the experimental Raman and infrared frequencies in Table 
I for reference. For NaCo02 with only one Na site occupied, symmetry analysis leads to 
the Raman-active modes T Raman = A\ g + E\ g + 2E 2g and the infrared-active optical phonon 
modes Tjpt = 2A 2u + 2E\ U . The experimental Raman frequencies listed in Table I generally 
agree with our calculated frequencies well. Among the four infrared-active modes, only the 
frequency of the hard E lu mode is reported in experiments. The frequency of Nao.57Co02 
reported by Lupi et al.— is near our calculated TO frequency of E lu mode. It is very strange 
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that infrared spectra by Wang et al.— gave four different infrared frequencies for doubly- 
degenerated in-plane Ei„ mode in metallic Nao.7Co02, where no LO-TO splitting exists. 
There should be only two frequencies even considering the partial Na occupation of the 2b 
and 2d site at the same time. We argue that part of these frequencies may be contributed by 
defects or grain surface, since the IR reflectivity data are very sensitive to surface treatment. 
In a word, the IR modes remain to be verified and understood. The peak of optical phonon 
at 161.3 cm" 1 (20 meV) in neutron inelastic scattering experiment is near the frequency 
of the soft E 2g mode (mode 1) here. 

IV. CONCLUSION 

In conclusion, we have calculated the lattice dynamics of NaCoC>2. Most of the phonon 
modes are only little affected by the Na site occupancy, and the calculated zone-center 
phonon frequencies generally agree with those by recent Raman and neutron scattering 
spectra on corresponding Na de- intercalated and/or hydrated compounds well. Therefore 
our calculations may also be helpful for understanding the lattice dynamics of the cobalt 
oxide superconductor. However, to clarify the phonon contribution in the superconductivity 
requires much further work on the lattice dynamics and electron-phonon interaction of the 
superconductor itself. Experimental data on infrared-active modes are difficult to compare 
with our calculated results now. In addition to phonon frequencies, other useful data such 
as Born effective charges and dielectric constants of NaCoC>2 are also presented. 
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TABLE I: Zone center optical phonon modes of NaCo02- The three groups separated by horizontal 
lines are Raman-active, infrared-active and silent modes respectively. The involved atoms in each 
mode and their vibration directions (/ / for in-plane, _L for out-of-plane) are listed. The atomic 
displacements in parentheses are obtained by dividing the normalized eigenvector components by 
the square root of the atomic mass, and with the unit of 10~ 2 for convenience of reading. These 
displacements are only refer to geometry A, since they are very similar to those of geometry B. 
f A and f B stand for calculated frequencies for geometries A and B respectively, and f EXP stands 
for experimental frequencies. For the four infrared-active modes, both the LO and TO frequencies 
are listed. 



Mode Symmetry 


Vibrations 


f A (cm~ 1 ) 


Acm" 1 ) 


/^(cm- 1 ) 


1 




Na//(0.34) O//(0.04) 


172.9 


185.7 






2 




O//(0.29) 


477.1 


482.1 


458 a . 


, 480 b , 469 c 


3 


"is 


Na//(0.05) O//(0.29) 


483.7 


489.8 


494 a 




4 


Ai„ 


O_l(0.29) 


608.0 


604.6 


574 a . 


,598 b ,582 c ,588 d 





TP, 


Na//(0.31) Co//(0.09) O//(0.05) 


201.2 


216.5 












206.3 


220.4 






6 


A 2u 


Na_L(0.31) Co_l(0.09) O_l(0.05) 


397.5 


337.0 












418.1 


361.6 






7 


A 2 « 


Na_L(0.03) Co_l(0.12) O_l(0.24) 


569.8 


566.5 












618.7 


615.1 






8 




Na//(0.03) Co//(0.12) O//(0.24) 


586.7 


590.0 


570 e , 


(505, 530, 560, 575) / 








638.5 


638.6 






9 




Co//(0.18) O//(0.16) 


88.0 


95.2 






10 


B2m 


Co ± (0.18) j. (0.15) 


197.0 


172.1 






11 


Big 


Na_L(0.33) O_l(0.10) 


351.7 


309.9 






12 


E>2u 


Co//(0.12) O//(0.24) 


582.8 


585.6 






13 


B2m 


Coj.(0.11) O ± (0.25) 


616.2 


610.0 






14 


Big 


Na ± (0.11) O_l(0.28) 


622.4 


616.0 







a Raman frequency of Nao.7Co02- 

6 Raman frequency^ of Na x Co0 2 - yH20 powder samples with x = 0.35. 
c Raman frequency^ of Nao.3CoO2-l.3H2O single crystals. 
d Raman frequency^ of Nao.7Co02 single crystals. 
e Infrared frequency 14 ; of Nao.57Co02- 



TABLE II: The Born effective charges of NaCoC>2 atoms. See text for the definition of Z//, Zj_ 
and Z. Superscripts A and B refer to data for geometries A and B. 





Z // 




Z 


Z A (Na) 0.87 


1.37 


1.03 


Z A (Co) 


2.49 


0.87 


1.95 


Z A {0) 


-1.68 


-1.12 


-1.49 


Z s (Na) 


0.85 


1.40 


1.04 


Z B (Co) 


2.41 


0.78 


1.86 


Z B (0) 


-1.63 


-1.09 


-1.45 
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FIG. 1: (color online) (a) Structure model of geometry A. The smallest balls represent cobalt 
atoms, and the balls in red bonding to Co are oxygen atoms. Between CoO layers, there exist 
sodium (big purple ball) layers, (b) Structure model of geometry B. (c) Electronic band structure 
and density of states of NaCo02 • 
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FIG. 2: Phonon band structure and density of states of NaCoC>2 for (a) geometry A and (b) 
geometry B respectively. 




